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ABSTRACT 

Using Monte Carlo simulations, we demonstrate photopion production from 
Fermi-accelerated protons and the resulting neutrino production in gamma-ray 
bursts. Unless internal shocks occur at quite large distance from the center, ul¬ 
tra high-energy protons are depleted by photopion production and synchrotron 
radiation. Internal shocks at hducial distance cause neutrino bursts, which ac¬ 
company gamma-ray bursts originating from electromagnetic cascades. 

Subject headings: gamma rays: bursts—cosmic rays—neutrinos 


1. INTRODUCTION 


The rapid time variabilities and the compactness problem (see, e.g., a review by Piran 
1999) suggest that gamma-ray bursts (GRBs) should arise from internal shocks within rel¬ 
ativistic flows. In the standard model a strong magnetic held is generated, and electrons 
are Fermi-accelerated in shocked regions. Because of the high Lorentz factor of the hows, 
synchrotron radiation by accelerated electrons is blue-shifted, and observed as gamma-ray 
photons. 

The physical conditions in the shocked region imply (Waxman 1995) that protons may 
be also Fermi-accelerated to energies > 10^° eV. If the total energy of protons, accelerated 
over the energy range 10^®-10^^ eV, is comparable to the energy of the gamma-rays, it can 
explain the hux of ultra-high-energy cosmic rays (UHECRs). This is an attractive idea for 
the cosmic-ray source. 

High energy protons in the GRB photon held can create high-energy neutrinos via 
photopion production (Waxman and Bahcall 1997, 1999; Guetta et al. 2001; Dermer and 
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Atoyan 2003; Guetta et al. 2004). If the above prediction is true, a Ikm^ neutrino detector 
may detect neutrinos correlated with GRBs in the near future (Waxman and Bahcall 1997; 
Wick et ah 2004). However, if the proton photopion “optical depth” is much higher than 
unity, protons lose their energies in GRB sources. In this case GRBs cannot be the sources of 
UHEGRs. Asano and Takahara (2003) (hereafter AT03) analytically show that high-energy 
protons cool rapidly via photopion production for fiducial parameters in GRB physics and 
that most of the proton energy converts to photons via 7r° decay, pion synchrotron, muon 
synchrotron, and the electromagnetic cascade of electron-positron pairs. 

In this paper, using the Monte Garlo method, we simulate numerically the proton cool¬ 
ing and resulting neutrino emission in order to conhrm quantitatively whether GRBs can 
be sources of UHEGRs or not. Our method pursues energy loss processes of each nucleon 
interacting with photons and the magnetic held during the dynamical timescale. Therefore, 
we can treat energy loss processes correctly even in the case that a nucleon creates pions 
multiple times in the dynamical timescale. Dermer and Atoyan (2003) carried out similar 
calculations using a different method from ours. Our purpose is to obtain a physical con¬ 
dition to produce UHEGRs, while Dermer and Atoyan (2003) mainly discussed neutrino 
production. In §2 we explain our model and method. The physical situation we consider is 
clearly shown. The numerical results are in §3. An analytically detailed description of the 
results can be found in AT03. Our conclusions are summarized in §4. 


2. MODEL AND METHOD 

We simulate nucleon cooling in the GRB photon held by the Monte Garlo method. We 
adopt a relatively low GRB energy and a small magnetic held for UHEGR production to be 
as advantageous as possible. The situation we consider is as follows. 

In the standard model gamma-ray photons are emitted from shocked shells moving 
with the bulk Lorentz factor T > 100 (see, e.g., a review by Piran 1999). In most cases the 
observed isotropic energies of GRBs are larger than 10^^ ergs (Frail et al. 2001), and the 
number of pulses per burst is less than 10 (Mitrofanov et al. 1998). According to the above 
facts, we consider one shell of T = 300, in which gamma-ray photons are produced, and the 
total energy of those photons is = 10^^ ergs in spherical symmetry. The shell width is 
constrained by the observed variability time-scale, 5t > 1 ms, and the duration time, ~ 10 s. 
Therefore, we assume that the shell width in the shell rest frame I is in the range of lO^'^-lO^"^ 
cm. This implies the average luminosity L = 9 x 10^®-9 x 10^^ ergs s“^. The photon energy 
density in the shell rest frame depends on the radius R ~ r‘^c6t where 

an internal shock occurs. In order to accelerate protons to energies > 10^° eV, the internal 
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shocks should occur at radii > 10^^ cm (Waxman 1995). Therefore, we consider cases with 
R > 10^^ cm only. The energy density of magnetic held Ub = B‘^/87r is assumed to be 

In our simulation, adopting the typical observed GRB spectra, the photon number 
density in the energy range + de^ in the shell rest frame is set at n{e^) oc for 1 eV 
< < 1 keV and for 1 keV < e.y < 10 MeV. The break energy 1 keV corresponds to 

300 keV in the observer frame. For e.y < 1 eV the synchrotron self absorption may be crucial 
(Granot et al. 2000), while the pair absorption may be crucial for > 10 MeV (e.g., see 
AT03). Thus, we neglect photons in such energy regions. 

Since GRBs are possible sources of UHEGRs with enegies > 10^° eV, we assume Fermi 
accelerated protons are injected in the energy range of ep = 10^°-10^® eV in the shell rest 
frame as iV(ep) oc Although the lowest energy of Fermi-accelerated particles is unknown, 
it is natural to consider that accelerated protons are distributed to quite a lower energy. 

Relativistic protons (cp = ypTUpC^ ^ rripC^) in the GRB photon held create pions and 
lose their energies via photopion production. In the shell rest frame the photon distribution 
is isotropic. Then, the time scale of photopion production, is written as 

^^^(7p) = 27rc j d/i(l-/3p/r) j (1) 

where (3^ = y^l — l/7p — 1 and /i is the cosine of the photon incident angle. From Stecker 
(1968), based on experiments, we approximate the photopion production cross section by 
a broken power-law prohle as (JTr{x) = 5 x 10“^®(x/590)^'^cm^ for 290 < y < 590 and 
<^ 7 r(x) = 5 X 10“^®(x/590)“°'^cm^ for 590 < x < 9800, where ymeC^ is the photon energy in 
the proton rest frame. The energy of a photon is expressed as = xmeC^/ 7 p(l — /3p/i). 

Secondary neutrons produced in the channel with isospin hip, py ^ nn~^, also create 
pions on almost the same time-scale as We use equation (1) for both protons and neutrons. 
In the A-approximation (Stecker 1973) the creation ratio : 7r° is 1 : 2. In the power- 
law target photon held, however, the ratio : 7r° is close to 2 : 1 (Rachen and Meszaros 
1998), because the charged pion production rate increases away from the A-resonance. So, 
we adopt = 2 : 1 in this simulation. For simplicity the nucleon inelasticity is hxed as 

K = Acp/cp = Acn/en = 0.2, and the process of multiple pion production, which enhances 
the nucleon cooling efficiency, is neglected. 

Extremely high energy protons can cool via synchrotron radiation with the power per 
unit photon energy 


dW 
dt dd'y 


sin a 
hmpC^ 


F{x), 


( 2 ) 
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where x = 2mpCe^/{‘i'y'^heB sin a) and F{x) is the synchrotron function (Rybicki and Light- 
man 1979). The pitch angle a is determined from random numbers in our Monte Carlo simu¬ 
lation. The cooling time scale due to synchrotron radiation is fsY = 3(7p—l)mpC®/2e^i?^7p/5p sin^ a. 

Another possible mechanism of the proton cooling is inverse Compton (IC) scattering. 

The scattering rate is obtained from the Klein-Nishina cross section ctkn (Rybicki and Light- 
man 1979) by exchanging mass nie rrip. In the proton rest frame, when an photon with 
energy Cq is scattered, the energy of the scattered photon is 


e 


/ 

1 


^ ~ T'-s) 


( 3 ) 


where /is is the cosine of the scattering angle. The probability distribution of the scattering 
angle is 


P oc — 
^0 




( 4 ) 


The Lorentz transformation yields the scattered photon energy ei in the shell rest frame. 
The cooling time-scale due to IC is 


^ic (bp) — 


n[e^ 


(7p - l)mpc‘ 


:27rc / (i/i(l -/3p/i) / iie^(ei - e^)^;^(TKN(eo) 




( 5 ) 


where Cq = 7 pe..y(l — /5p/i). To obtain tic, the numerical integration is carried out with /is 
determined from equation (4). This numerical integration process also gives the IC spectra 
at once. 


We pursue 2.7 x 10^ nucleons in each computation. Given a nucleon energy, the cooling 
time scales F, tsY, and tic are determined. Let us denote the cooling time scale due to 
photon emission by ty, which satishes = tgy -l- Assuming that nucleons interact 
with photons and the magnetic held during the dynamical time scale = l/c, we pursue 
the nucleon cooling process by a time step At = min(t 7 r, t^)/10 and p/lO for protons and 
neutrons, respectively. The time step will change as nucleons lose their energies. The pion 
production probability in each time step is 1 — exp(—At/t^). If nucleons create charged 
pions, the species of the nucleons is changed {p n). After td s have passed, nucleons are 
considered to escape from the shell. 

The neutron-decay emissions will produce nonthermal emission signatures. However, 
we can distinguish those emissions from the prompt burst, since the hux of neutrinos and 
photons originated from the neutron-decay is faint because of the long decay time-scale. We 
do not discuss those emissions below. 
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3. RESULTS 
3.1. Proton cooling 

The time-scales tsy, and tic are proportional to R^l, while the dynamical time-scale 
td oc I (see AT03). The proton cooling process is determined by the ratios = td/t^, 
Tsy = td/tsY, and Tie = td/tic, which are proportional to R~^. Thus, the proton cooling 
process is independent of 1. Under the situation we consider, the proton cooling efficiency 
depends on the total photon energy and the radius R, rather than on the luminosity 
L oc E^/l. 

In Figure 1 we plot the ratios Tt,,tsy,tic for R = 10^^ cm and E^ = 10^^ erg. From 
the condition to create pions > 0.2GeV^, the behavior of is explained as follows. For 
Fep < 6 X 10^® eV, protons produce pions with photons whose energies are above the break 
energy Fe^ = 300 keV [n{e^) oc In this case the number of the interacting photons 

is proportional to ~ Cp^, which implies oc Cp^. Above 6 x 10^® eV protons interact with 
photons in a low-energy region [n(e.y) oc As a result, is nearly constant in this 

energy region. It is easily conhrmed that the simple approximation in Waxman and Bahcall 
(1997), which was adopted in simulations in Guetta et ah (2001, 2004), produces a close 
value to our numerical estimate, r,r- Since Waxman and Bahcall (1997) adopted a very low 
energy, E^ = L6t = 10"^® ergs for a shell, the pion production efficiency in their case becomes 
lower than ours. Beyond 6 x 10^® eV, the threshold energy of target photons to create pions 
is below the cut-off energy (0.3 keV in the observer frame). As a result, decreases above 
this energy. 

For Fcp > 3 X 10^° eV, the proton synchrotron cooling dominates photopion production. 
Although the photon energy density is larger than the magnetic energy density {Ub = 0.177,,), 
IG emission is negligible because of the Klein-Nishina suppression. 

Figure 1 indicates that all protons above 10^^ eV may lose their energies for R = 10^^ 
cm. The result of the numerical simulation shown in Figure 2 is consistent with the above 
discussion. Nucleons with energies > 10^® eV are extinguished. The energy protons lose 
is mainly transferred to pions and synchrotron photons (see Table 1). The highest energy 
protons (> 6 X 10^® eV) cool via synchrotron radiation, so that the pion distribution has a 
break at ~ 10^® eV. The low-energy tail of the synchrotron spectrum exhibits n{e.y) oc as 
is usually seen in the fast cooling case (Granot et ah 2000). The weak peak at ~ 10^ eV in the 
synchrotron spectrum may be due to low-energy protons whose photopion “optical depth” 
is small enough. The IG spectrum is complex because of the Klein-Nishina suppression. 

Even for R = 10^^'® cm (Fig. 3), nucleons > 10^® eV are depleted by 2 orders of 
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magnitude. Figure 4 shows the case for R = 10^"^ cm. In this case GRBs can be sources 
of UHECRs, although nucleons > eV are reduced by the synchrotron cooling. The 

maximum value of is ~ 1 because of the low photon density. The highest energy protons 
create a few pions and lose a few tens of percent of their energies. 


3.2. Neutrinos 

Using the same method, we simulate behavior of charged pions and muons until they 
decay into positrons (electrons) and neutrinos. The life times of pions and muons, = 
77 rTjr,o and (T,r,o = 2.6 x 10“® s and = 2.2 x 10“® s are the life times of 

pions and muons at rest, respectively) are independent of the shell width 1. Therefore, the 
energies charged pions and muons lose before their decays depend on I differently from the 
proton cooling. 

When a pion decays 7r+(7r“) ^ /i+z/^(/i“z7^), the energy fraction a muon obtains is 
~ ~ 0.76. The rest of the energy goes to a neutrino. We assume that the energy 

of a muon at its decay will be shared equally by a positron (electron), a neutrino, and an 
anti-neutrino. We neglect neutrinos due to neutron decay n pe~i7e, whose time-scale is 
much longer than t^- 

Figure 5 shows the resulting spectra of neutrinos and photons emitted by pions and 
muons for R = 10^^ cm and I = 10^° cm. Pions and muons above the energy where the life 
time is comparable to the cooling time t.y or dynamical time td (Waxman and Bahcall 1997; 
AT03) lose their energy before their decay. For R = 10^^ cm and I = 10^'^ cm, one half of 
the initial energy of the pions is expended via synchrotron radiation, and muons also lose 
their energy via synchrotron (~ 35%) and IC (~ 17%) radiation. Neutrinos are distributed 
as ny{ey) oc with a low- and high-energy cut-off. 

The high-energy cut-off is determined by the life-time and the cooling time of pions 
(muons). Since a larger I leads to a lower Ub, the high-energy cut-off is roughly proportional 
to (see Fig. 6). Wider shell decreases the contribution of IC photons in comparison with 
synchrotron radiation, because of the Klein-Nishina suppression. 

The low-energy cut-off is determined by the minimum energy of pions produced from 
protons, which depends on only R, as was discussed in §3.1. Both the low-energy and high- 
energy cut-off increase roughly oc R (see Fig. 7). If we observe neutrinos from GRBs, the 
spectra give us parameters such as I and R independently of gamma-ray observation. 

Figure 8 shows the ratio of the neutrino energy to the initial proton energy (> 10^'^ eV 
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at injection). For R < 10^^ cm, ~ 10% of the proton energy will be converted to a neutrino 
burst. If we consider only muon neutrinos, which will be detected with a kilometer-scale 
detector, their energy fraction may be a few percent of the accelerated protons. This result 
is consistent with the prediction of Dermer and Atoyan (2003). 

The neutrino energy is a part of the energy protons lose via photopion production and 
synchrotron radiation. The “rest energy” (the energies of neutral pions, photons, electrons, 
and positrons) will be converted to lower energy photons (< 10® eV) via electromagnetic 
cascades (see AT03), because higher energy photons (> 10® eV) will be absorbed by electron- 
positron pair creation. The total energy of such photons is quite larger than the neutrino 
energy (see Fig. 9). If the energy of Fermi-accelerated protons above 3 x 10^® eV is compa¬ 
rable to the GRB photon energy the energy of protons > 10^® eV is 4.5E^. Thus, the 
energy of protons may be 1-5 E.y. In any case the energy of nucleons above 3 x 10^® eV will 
be lost for R < 10^^ cm. As long as the energy of UHECRs is assumed to be comparable to 
the GRB energy, gamma-ray photons originating from protons may be a large contiribution 
to the GRB photons. 


4. CONCLUSIONS AND DISCUSSION 

Our simulations have shown that internal shocks should occur at radii > 10^^ cm to 
generate UHECRs, although siginihcant GRB pulses have shorter time-scales of < 1 s (Norris 
et ah 1996). On the other hand, neutrino bursts are expected for R < 10^"^ cm, while ultra- 
high-energy protons are depleted. The energy fraction of the neutrinos and neutrino spectra 
have been calculated for various R and /. Neutrino spectra give us information on the physical 
condition of GRBs, although only the most powerful bursts, which are brighter than 10^® 
erg or occur at x <~ 0.1, produce detectable neutrino bursts with a kilometer-scale detector 
(Dermer and Atoyan 2003). 

Smaller E^ and larger R are favorable to generate UHECRs. At i? = 10^^ cm, the 
most probable shell width is / ~ R/T ~ 10^^ cm. In this case the magnetic held becomes 
8 X 10® G assuming E^ = 10®^ ergs and Ub = 0.17%. This strength seems rather weak 
to produce observed gamma-ray photons. Since there are some unsolved problems and 
ambigious points in the theory of GRB spectra (Preece et ah 1998; Totani 1999; Asano and 
Kobayashi 2003), we have simply adopted a typical observed spectra without considering 
any physical conditions. It is challeging to hnd a condition to make both UHECRs and 
typical gamma-ray spectra be consistent. 

To produce both UHECRs and a neutrino burst at the same time, the radius is limitted 



to around R ~ 10^^ cm in our case. Throughout this paper, we have adopted a relatively 
small GRB energy = 10^^ ergs as an optimistic case. For brighter bursts {E^ > 10^^ ergs), 
stronger magnetic helds {Ub ^ U^), and smaller Lorentz factors (F < 300), internal shocks 
should occur at more distant radii to generate UHECRs. The threshold radius behaves as 
roughly oc under our assumption. 

If protons are accelerated to very high energies at both the internal shock phase and 
afterglow phase, both UHECRs and neutrino burst may be produced. In this case the 
neutrino bursts may be originated from internal shocks, while most of the UHECRs are 
produced at the afterglow phase. 

Throughout this paper the internal shock model has been applied. However, the external 
shock model may be applicable for smooth prohle CRBs. According to Dermer (2002), the 
neutrino flux from individual CRBs is far too weak in the external shock model, while 
neutrino bursts occur in the internal shock model as shown in this paper. In the external 
shock model, ultra-high-energy neutrons are produced and travel ~ 100 kpc before their 
numbers are depleted by /5-decay. Therefore, the external shock model predicts detection 
of neutron /5-decay halo emission at optical and radio frequencies with luminosity ~ 10^^ 
ergs s“^. On the other hand UHECRs may be extinguished for bright CRBs in the internal 
shock model so that neutron /5-decay halos are not generated. Euture neutrino, optical, and 
radio observations may provide tests of the two scenarios. Measurement of ion composition 
of UHECRs from supernova remnants with instruments such as Auger may be helpful for 
investigating the origin of UHECRs. 

Neutrino bursts surely accompany gamma-ray photons originating from pion cascades. 
Those photons may be not negligible in comparison with the “original CRB photons.” In 
this case, as AT03 discussed, CRB photons may arise from pion cascades, rather than from 
initially existing electrons. 


The author thanks the referee for useful comments. 
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Fig. 1.— Proton cooling efficiencies r.^-, t'sy, Pc for R = 10^^ cm and = 10^^ ergs. The 
proton energies are measured in the observer frame. 
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Fig. 2.— Spectra of GRB photons, protons at injection, created pions, photons emitted by 
protons {dashed lines)^ and cooled nncleons (protons and nentrons) for R = 10^^ cm. 







13 


Log[£^n(E)] 
(eV /crn^) 



Fig. 3.— Same as Fig. 2, but for R = cm. 
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Fig. 4.— Same as Fig. 2, but for R = 10^^ cm. 
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Fig. 5.— Spectra of neutrinos {solid lines) from pion decay and and muon decay 
{ue, Fe, and and synchrotron and IC photons emitted by pions {thin dashed lines) and 
muons {thick dashed lines) for R = 10^^ cm and I = 10^*^ cm. 
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Fig. 6.— Same as Fig. 5, but for I = 10^^ cm. 
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Fig. 7.— Spectra of neutrinos from pion decay for I = 10^° cm. The radius R changes from 
10^3 to cm. 
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Fig. 8.— Ration of neutrino energy to the initial proton energy. Dashed lines are for 
neutrinos originated from pions only, and solid lines are the total fraction. The shell width is 
assumed to be 10^° {open circles), 10^^ {rectangles), and 10^^ {filled circles) cm, respectively. 
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Fig. 9.— Ratio of “rest energy” (see text), which will be converted to gamma-ray photons, 
to the initial proton energy. The shell width is assumed to be 10^° {open circles), 10^^ 
(rectangles), and 10^^ (filled circles) cm, respectively. 
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Table 1: Energy Fraction Each Species of Particle Contributes to the Total Energy of Protons 
at Injection. 


R (cm) 

N 

71 

7 

10^3 

0.33 

0.57 

0.10 

1013.5 

0.43 

0.47 

0.10 

10'" 

0.73 

0.19 

0.08 

1014.5 

0.94 

0.03 

0.03 

lO's 

0.99 

0.002 

0.01 


Note. — The notations N, tt, and 7 represent cooled nucleons, created 7 r° and tt^, and photons emitted 
directly by protons, respectively. 



